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Abstract. We have extended our previous mean-field galactic dynamo model which included algebraic and dynamic 
alpha nonlinearities (Kleeorin et al. 2002), to include also a quenching of turbulent diffusivity. We readily obtain 
equilibrium states for the large-scale magnetic field in the local disc dynamo model, and these fields have strengths 
that are comparable to the equipartition field strength. We find that the algebraic nonlinearity alone (i.e. quenching 
of both the a effect and turbulent magnetic diffusion) cannot saturate the growth of the mean magnetic field; 
only the combined effect of algebraic and dynamic nonlinearities can limit the growth of the mean magnetic 
field. However, in contrast to our earlier work without quenching of the turbulent diffusivity, we cannot now find 
satisfactory solutions in the no-z approximation to the axisymmetric galactic dynamo problem. 
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1. Introduction 

Spiral galaxies possess large-scale magnetic fields whose 
spatial scale is comparable with galactic radii (see for re- 
view Beck et al., 1996). Galactic magnetic fields mainly 
lie in the galactic plane and the corresponding magnetic 
lines are usually roughly of spiral form. This form can be 
. substantially distorted in the presence of strong noncircu- 
' lar motions, e.g. in barred galaxies, see Beck et al. (1999, 
2002), Moss et al.(2001). 

Galactic magnetic fields are believed to originate in a 
galactic dynamo, driven by the joint action of the mean 
hydrodynamic helicity of interstellar turbulence and dif- 
ferential rotation. The linear stage of galactic dynamo ac- 
tion seems now to be well- understood, see, e.g., Ruzmaikin 
et al. (1988). The conventional approach to the nonlinear 
stage of the galactic dynamo is based on comparison of 
the relative intensity of three quantities participating in 
dynamo action, namely the differential rotation 6fl, turbu- 
lent diffusivity rjx and a-effect, each of which can be asso- 
ciated with a typical velocity: 200 km s^^, 10 km s^^ and 
1 km respectively. Because the typical velocity associ- 
ated with a is the smallest, the mean hydrodynamic helic- 
ity is believed to be the weakest part of the self-excitation 



Send offprint requests to: I. Rogachevskii 



chain, and a scenario of nonlinear dynamo limitation via 
a-quenching is the most commonly adopted. 

A simple version of a-quenching prescribes the mean 
hydrodynamic helicity to be a decreasing function of mean 
magnetic field strength B. The critical magnetic field 
strength i3cr, at which quenching becomes significant, 
is estimated conventionally from equipartition with the 
kinetic energy of interstellar turbulence, B^^ ^ Airpv^. 
When applied to specific galaxies, this picture results in 
robust magnetic field models which are compatible with 
observations. However the picture is obviously oversim- 
plified and various attempts to suggest a more developed 
version of nonlinear galactic dynamo theory have been un- 
dertaken by several authors. In particular, Vainstein and 
Cattaneo (1992) and Gruzinov and Diamond (1995) have 
claimed that in fact _Bcr is much lower then the equiparti- 
tion value, e.g. Bcr = Scqi?m~^/^, where Rm is the mag- 
netic Reynolds number. In galaxies, magnetic Reynolds 
numbers are very large, Rm « 10^ even if the ambipolar 
diffusivity coefficient is used, so it was claimed that dy- 
namo action saturates at a magnetic field strength that 
is much lower than both the equipartition value, and also 
the large-scale field strengths observed in nearby spiral 
galaxies. This result follows from investigations that ac- 
cept the idea of magnetic helicity conservation. The galac- 
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tic dynamo produces a large-scale magnetic field with 
nonvanishing magnetic helicity, and when considering de- 
tailed magnetic helicity conservation in a given volume, 
the above 'catastrophic' estimate for B„ results. 

The evolution of magnetic helicity appears however to 
be a more complicated process than can simply be de- 
scribed by a balance of magnetic helicity in a given vol- 
ume, and it is necessary to take into account, as for the 
evolution of other conserved quantities, transport by the 
fluid flows including turbulent transport of magnetic helic- 
ity through the galactic boundaries and the destruction of 
mean magnetic helicity by turbulent diffusion. The gov- 
erning equation for magnetic helicity was suggested by 
Kleeorin & Ruzmaikin (1982; see also the discussion by 
Zeldovich et al. 1983) for an isotropic turbulence, and in- 
vestigated by Kleeorin et al. (1995) for stellar dynamos, 
and self-consistently derived by Kleeorin & Rogachevskii 
(1999) for an arbitrary anisotropic turbulence. A quan- 
titative model for the flux of magnetic helicity was pro- 
posed by Kleeorin & Rogachevskii (1999) and Kleeorin et 
al. (2000). Note that Schmalz & Stix (1991), Covas et al. 
(1998) and Blackman & Brandenburg (2002) have also in- 
vestigated related solar dynamo models that included a 
dynamical equation describing the evolution of magnetic 
helicity. Magnetic helicity transport through the bound- 
ary of a dynamo region is reported by Chae (2001) to be 
observable at the solar surface. The role of a flux of mag- 
netic helicity in the dynamics of the mean magnetic field 
in accretion discs was also discussed by Vishniac & Cho 
(2001). 

The equation governing the magnetic helicity is much 
more complex than the conventional parametrization used 
to represent a-quenching. Kleeorin et al. (2000) suggested 
that a nonlinear galactic dynamo governed by a consis- 
tently derived equation for magnetic helicity results in a 
steady magnetic field comparable with the cquipartition 
magnetic field estimate, and Kleeorin et al. (2002) demon- 
strated that a detailed galactic dynamo model based on 
the equation under discussion gives rc^sults very similar to 
one based on conventional a-qucncliing. In other words, 
the real physical description of the nonlinear stage of 
galactic dynamo is quite complicated but, if we are inter- 
ested in pragmatic results only, an adequate description 
can be found by using only a conventional a-quenching. 

We stress that the scenario of Kleeorin et al. (2000, 
2002) does not include all possible types of nonlinear pro- 
cesses which can occur at the nonlinear dynamo stage (see, 
e.g., Brandenburg & Subramanian 2000; Brandenburg & 
Dobler 2001), but rather is restricted by a minimal number 
of processes involved in the magnetic helicity conservation. 
In particular, we consider transport of mean helicity of 
small-scale magnetic field only and note that this helicity 
may be transported out of the galactic disc without sig- 
nificant losses of large-scale magnetic field. Moreover, this 
picture is formulated in terms of mean-field electrodynam- 
ics, which gives a natural constraint on the description of 
small-scale details. When attention is focussed on these 
details (e.g. Brandenburg & Sokoloff 2002; Blackman & 



Brandenburg 2002), the mean-field description should be 
considered as a parametrization of the turbulence. 

In the spirit of the basic ideas about the nonlinear 
saturation of galactic dynamos, the analysis presented by 
Kleeorin et al. (2000, 2002) was restricted to the evolution 
of a only, while a detailed simulation (e.g. Brandenburg 
& Sokoloff 2002) also demonstrates a quenching of the 
turbulent magnetic diffusivity. A quantitative model for 
a nonlinear quenching of turbulent magnetic diffusivity 
has been recently suggested by Rogachevskii & Kleeorin 
(2001). 

The aim of the present paper is to include a turbu- 
lent magnetic diffusivity quenching into the mean-field 
dynamo equations. This effect is expected to be quite 
modest. Speaking pragmatically, we do not know the tur- 
bulent magnetic diffusivity of interstellar turbulence well 
enough to recognize its saturation, by some dozens of per- 
cent. However, our analysis below demonstrates that the 
problem is not restricted by some specific variation of 
magnetic diffusivity coefficient. Because of nonlinear ef- 
fects, the turbulent magnetic diffusion coefficients for the 
two basic magnetic field components, i.e. poloidal and 
toroidal, become different (see Rogachevskii & Kleeorin 
2001). Since spatial magnetic field distribution is nonuni- 
form, a nonuniform magnetic diffusivity saturation arises, 
that results in new terms in the dynamo equations govern- 
ing the nonlinear magnetic field evolution. In general, the 
situation appears to be less trivial then might be thought 
initially, and a quantitative analysis of a specific dynamo 
model becomes desirable. Below we present results of the 
corresponding analysis and numerical simulations. 

2. The mean-field equations 

The mean-field dynamo equation (e.g. Krause & Radler 
1980) is 

— =Vx{VxB + e-r]VxB), (1) 

where V is a mean velocity {e.g., the differential rota- 
tion), r] is the magnetic diffusion due to the electrical 
conductivity of fluid together with ambipolar diffusion, 
6 — (u X b) is the mean electromotive force, u and b 
are fluctuations of the velocity and magnetic field, respec- 
tively, angular brackets denote averaging over an ensem- 
ble of fluctuations. When a small-scale background turbu- 
lence (i.e. turbulence with a zero mean magnetic field) is 
isotropic and the energy of small-scale magnetic fluctua- 
tions of the background turbulence is much smaller than 
that of the kinetic energy, the mean electromotive force is 
given by 

£i = aij {B)Bj + ( V^fi^(B) x B)i - mj (B) (V x B)j , (2) 
where 

Vij (B) = [r]A (B) - TIP (B)] + rjp (B) Bij , (3) 
V^^iB) = lvv{B)^ , (4) 



N. Kleeorin et al.: Nonlinear magnetic diffusion and magnetic helicity transport in galactic dynamos 



3 



(see Rogachevskii & Kleeorin, 2001), Bij = BiBj/B^, 
a^j{B) = iiv{B) = t^a{B)-t^b{B)-i^p[B) and 

a{B), ?7^(B), i]b{B), 'qp{B) are determined by Eqs. (^, 
(|ll|), @ and (0), respectively. 

2.1. The local thin-disc dynamo problem 

We begin by considering the simplest local disc dynamo 
problem. Using cylindrical polar coordinates r,(j), z, from 
Eqs. ([|)-(|4|) we obtain the equations for the mean radial 
field Br — Rabr and toroidal field B^ for the local thin- 
disc axisymmetric cki7-dynamo problem as 

^ ^ -HB)B^y + iVA{BWJ - iVA{B)bry , (5) 
dB. 



dt 



DK + {MB)B')' 



(6) 



(Ruzmaikin et al. 1988, Rogachevskii & Kleeorin 2001). 
Here a prime denotes d/dz^ Q;(B) is the total nonlin- 
ear a effect, r]A{B) and rjB{B) are the nonlinear turbu- 
lent magnetic diffusion coefficients of poloidal and toroidal 
mean magnetic fields, and the nonlinear function Va{B) = 
[r]A[B) - riB{B)][\nB)' , with B = \B\. These nonlineari- 
ties are specified and discussed in the next Section. 

We adopt here the standard dimensionless form of the 
galactic dynamo equations from Ruzmaikin et al. (1988); 
in particular, length is measured in units of the disc thick- 
ness h, time in units of h? /rjT and B is measured in units 
of the equipartition energy i?cq = \/47r/9u», a is measured 
in units of a* (the maximum value of the hydrodynamic 
part of the a effect), the nonlinear turbulent magnetic dif- 
fusion coefficients riA{B) and riB{B) are measured in units 
of i]T. We define Ra — ha^,/r]T, Rui — r {dQ / dr) h'^ / rjT , 
and the dynamo number D = R^jRa, where / is the max- 
imum scale of the turbulent motions, Rm = lUi,/r] is the 
magnetic Reynolds number. Also is the characteristic 
turbulent velocity at the scale I, p is the gas density, and 
the characteristic value of the turbulent magnetic diffusiv- 
ity riT = 

Throughout this paper, unless otherwise stated, we 
measure magnetic field in units of the equipartition value. 
Here we assumed that the background turbulence (i.e., the 
turbulence with a zero mean magnetic field) is isotropic 
and has only velocity fluctuations, even though a nonzero 
mean magnetic field can be expected to produce an 
anisotropy of turbulence from the generated magnetic fluc- 
tuations. 



3. The nonlinearities 

3.1. The nonlinear a effect 

The total nonlinear a effect is given by 
a{B) = a" a"* , 



(7) 



where is the hydrodynamic part of the a effect, 
and a™ is the magnetic part of the a effect. These 



quantities are determined by the corresponding helici- 
ties and quenching functions, (jtv^B) and 4>m{B). In par- 
ticular, a'' = x>.(B), a™ = x%B)(t>„,{B) and = 
-(r/3)( u ■ (Vxm) ), x" = (T/127rp)(6 • (Vx6)) is re- 
lated with current helicity, where r is the correlation time 
of turbulent velocity field and {u ■ (Vxi6)) is the hydro- 
dynamic helicity. Thus, 



a{B)^x''4>v{B)+x'{B)4>UB) 



(8) 



where the quenching functions (pviB) and (f>miB) are given 

by 



MB) = (l/7)[4(/.„,(B) + 3L(\/8B)] , 
M{B) = (3//32)(l-arctan(/3)//3) 



(9) 
(10) 



(see Rogachevskii & Kleeorin, 2000), where (3 — VSB and 
L{P) = 1-2(3^ + 2/3Mn(l -I- f3~^). Thus 0^,(5) = 2/^"^ 
and <j),n{B) = 3//?2 for /3 > 1; and 0^,(5) = 1 - {6/ 5) (3^ 
and (j)rn{B) = 1 - (3/5)/32 for /3 < 1. The function x^iB) 
entering the magnetic part of the a effect is determined by 
the dynamical equation (p^. Here x^ &nd are measured 
in units of a*. 

The function 4iv{B) describes conventional quench- 
ing of the a effect. A simple form of such a quenching, 

= 1/(1 + B^), was introduced long ago (see, e.g., 
Iroshnikov, 1970). The splitting of the total a effect into 
the hydrodynamic {a") and magnetic (a™) parts was first 
suggested by Frisch et al. (1975). The magnetic part a™ 
includes two types of nonlinearity: the algebraic quench- 
ing described by the function (p„i{B) (see e.g. Field et al. 
1999; Rogachevskii & Kleeorin 2000, 2001) and the dy- 
namic nonlinearity which is determined by Eq. (p^). 

3.2. Nonlinear turbulent magnetic diffusion coefficients 
of the toroidal and poloidal mean magnetic fields 

The nonlinear turbulent magnetic diffusion coefficients of 
poloidal and toroidal mean magnetic fields riA{B) and 
riB{B), and the nonlinear function Va{B) = [riA{B) ~ 
riB{B)]{\ii B)' are given in dimensionless form by 

r,A{B) = A^{V2(3) + (l/2)A2(/3) , (11) 
r^BiB) = AiiV2p) - {l/2)A2{(3) + A2iV2(3) 

+ {f]/V2)'f{V2f3) , (12) 
Va{B) = [A2iP) ~ A2{V2I3) 

-(/3/^/2)vI/(V2/3)](lnB)', (13) 

(see Rogachevskii & Kleeorin, 2001), where '^{x) — 
{d/dx)[Ai{x) + (1/2)^2(2;)], and 

r,piB) = (l/2)A2(/3) . (14) 
The functions Ai{f3) and A2{P) are given by 
arctan f3 



AiW 
A2W 



f3 

arctan /3 

f3 



1 



5 \ 1 r/o^ 5 



1 



15 



15 
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When /3 ^ 1 these functions are given approximately by 
A,{P) = 1 - (2/5)/3^ , A^{P) = -(4/5)/32 , 
and for /3 S> 1 they are given by 

Ai{l3) = 37r/5/3 - 4/5/3^ , A2{I3) ^ -37r/5/3 + 14/5/3^ . 

The asymptotic formulas for the hmctions rjA{B), 
rjBiB) and Va(B) for /3 < 1 are 

rjA{B) = 1 - (6/5)/32 , T^siB) = 1 - (IS/S)/?^ , 
VaIB) - (12/5)/32(lnB)' , 

and for /3 > 1 



10/3 
20/3 



(2 



-1), 



The quenching of the a effect and the turbulent mag- 
netic diffusion are caused by the direct and indirect mod- 
ification of the electromotive force by the mean magnetic 
field. The indirect modification of the electromotive force 
is caused by the effect of the mean magnetic field on the 
velocity fluctuations (described by the tensors {muj)) and 
on the magnetic fluctuations (determined by the tensor 
(bibj)), while the direct modification is due to the effect of 
the mean magnetic field on the cross-helicity tensor {uibj) 
(see, e.g., Rogachevskii & Kleeorin 2000, 2001). 

3.3. The dynamical equation for the function 

The function x'^(B) entering the magnetic part of the a 
effect is determined by the dynamical equation 



(15) 



where F = Cx" (t)v{B)B'^T^A{B){V p) / p is the nonad- 
vective flux of the magnetic helicity which serves as 
an additional nonlinear source in the equation for x'^, 
Vx'^ is the advective flux of the magnetic helicity and 
— kVx'^ is the diffusive flux of the magnetic helicity 
(see Kleeorin & Rogachevskii, 1999; Kleeorin et al. 2000; 
2002), V — e^Ur is the differential rotation, and T — 
{l/3){l/h)'^Rm. Eq. (|l5|) was obtained using arguments 
based on the magnetic helicity conservation law (see 
Kleeorin & Rogachevskii, 1999). The function x'^ is pro- 
portional to the magnetic helicity, x'^ = X™/ {^^'^Vtp) (see 
e.g. Kleeorin & Rogachevskii, 1999), where x™ = {(J- ' b) 
is the magnetic helicity and a is the vector potential of 
small-scale magnetic field. Here we assume that the he- 
lical part of the vector potential a is a locally isotropic 
and homogeneous random field, which is a natural gauge 
condition used in our approach. Thus, Eq. ( p^ ) describes 



the evolution of magnetic helicity, i.e. its production, dis- 
sipation and transport. 

The turbulent diffusion of the magnetic helicity k plays 
an important role and can be interpreted as follows. The 
random flows existing in the interstellar medium consist 
of a combination of small-scale motions, which are af- 
fected by magnetic forces resulting in a steady-state of 
the dynamo, and a microturbulence which is supported 
by a strong random driver (supernovae explosions) which 
can be considered as independent of the galactic magnetic 
field. The large-scale magnetic field is smoothed over both 
kinds of turbulent fluctuations, while the small-scale mag- 
netic field is smoothed over microturbulent fluctuations 
only. It is the smoothing over the microturbulent fluctua- 
tions that gives the coefficient k. 

For galaxies the relaxation term jT is very small 
and can be dropped in spite of the fact the small yet fi- 
nite magnetic diffusion is required for the reconnection of 
magnetic field lines. For example, we will show below that 
the magnetic Reynolds number, i?m, does not enter into 
the steady state solution of Eq. ( p^ ) in the limit of very 
large i?m, because of the effect of the magnetic helicity 
flux. In particular, keeping the term x^ jT we obtain from 
Eq. (|l5|) that in a steady state 

_ a^+T, j?m(/>„(B)[.B- (Vx.B)-divFt] , . 

l + T,Rm(f>^{B)B^ ' ^ ^ 

where — l/u^, and Ft — F + Vx'^ ~ <^^X^ is the total 
flux of the magnetic helicity. In the limit of very large i?m, 
Eq. (|l6|) reads 



B ■ {VxB)-diYFt 



(17) 



This implies that in this limit the total a is independent 
of the magnetic Reynolds number. 

In the local approximation Eq. (|l^) reads: 

= 4(^y)'[rM(S)(i?^fe;-i?;6.)-«(S)B2 

+C^{\x^z)\UB)B'va{B))] + {Kixry ■ (18) 
oz 

Here we do not take into account any inhomogeneity of 
the turbulent magnetic diffusion at B = 0. The turbulent 
magnetic diffusion is inhomogencous due to inhomogene- 
ity of the mean magnetic field B. 

The turbulent magnetic diffusion k of the magnetic he- 
licity (and the function x*^) can depend on the mean mag- 
netic field. The nonlinear quenching of the turbulent mag- 
netic diffusion of the magnetic helicity is given by ^^(-B), 



l + Ai{V2f3) + ^A2{V2p) 



(19) 



For P <^ 1 we have approximately 
0.(5) = 1-^/32, 
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and for P ^ 1 we have 



The turbulent magnetic diffusion k{B) of the magnetic 
helicity is determined by the tensor T{uiUj). Thus, Eqs 
(0), (O)-(O) and (pq) contain the main nonlinearities. 



4. Equilibrium states of the local dynamo model 

4.1. Asymptotic expansions and an equilibrium solution 

We now present asymptotic expansions for a galactic dy- 
namo model determined by Eqs. (||), (||) and (^8|). For 
the afi-dynamo B « B^. This assumption is justified if 
\D\ ^ i?c(, i.e. \Ru,\ ^ 1. In a steady-state for fields of 
even parity with respect to the disc plane, Eqs. (^), (||) 
and (|l|) with K — gives 

[ilB{B)Bf + 2CDMB)VA{B)B^\x''iz)\ = . (20) 
The solution of Eq. (Um for negative D is given by 



/ G{B)dB = y/2C\D\ I y/Wmd~z 

Jo J\z\ 



(21) 



where G{B) = rjBiB)/[(j)^{B)TjA{B)B^]^/^ . For an arbi- 
trary profile |x''(z)|, negative dynamo number D and for 
B > l/VS, there is an explicit steady solution of this 
equation with the boundary conditions B^{z = 1) = 
and B'{z = 0) = 0, 



Biz) « i2/5)C\D\ 



V\rW\ds 



(22) 



where z is measured in the units of h, and we have used 
that for B > l/^/8 we have rjAiB) ~ 2/5/3, rjBiB) ~ 
3/5/3, (l3v{B) ~ 2//32. Here /3 = VSB. In a steady state 
br{z) — {riB{B)B'y /\D\. For the specific choice of the pro- 
file Ix^iz)] = sin2(7rz/2) we obtain 



B^ « —C\D\B,qCos' 



Br ~ — , „ , Bcq COS 



Ru 



nz 

Y 

■nz 

Y 



(23) 
(24) 



where we have now restored the dimensional factor B^q. 
The boundary conditions for Br are Br{z = 1) = and 
Br{z = 0) = 0. Note, however, that our asymptotic anal- 
ysis performed for B ^ l/VS is not valid in the vicinity 
of the point z = 1 because B{z = 1) = 0. 



4.2. Numerical solutions for the one-dimensional model 

We found solutions of Eqs. (||), (||) and (^ by step by 
step integration, from arbitrarily chosen initial conditions. 
Various properties of the solutions for the mean magnetic 
field and the z-profiles of the main nonlinearities - the 
a-effect, turbulent magnetic diffusion of the toroidal and 
poloidal fields - are illustrated in Figs. 1 - 6. It can be 




0.2 0.4 0.6 0.8 z 




Fig. 1. The z-dependence of solutions for the local model; 
h/l = 5; C = -0.1, K = 0.3 and |x"(2)| = sin2(7rz/2). 
The various curves indicate results with dynamo numbers 
D = -20 (dashed-dotted), D = -40 (dashed) and D = 
— 100 (solid), a) the toroidal magnetic field B^{z), b) the 
poloidal magnetic field Br{z), c) the total a-effect, a{z). 



seen clearly that the field strength (~ l-B^j) is typically of 
order 1 (equipartition) , and increases with \D\. 

New features were found in comparison with the re- 
sults of Kleeorin et al. (2002), where the quenching of the 
turbulent magnetic diffusion was not taken into account. 
In particular, the maximum of the toroidal magnetic field 
for \D\ > 40 is not located at z = but is shifted to 
z sa 0.5/i (see Figs, la and 3a). The equipartition toroidal 
field is attained for smaller values of the dynamo numbers 
and parameter C than in the model studied by Kleeorin 
et al (2002). The reason is that the asymptotic behav- 
ior of the steady-state solution is different in these two 



cases: in the present study B^ (x C \D\ B^q (see Eq. (p3[)) 
whereas when the quenching of the turbulent magnetic 
diffusion vanishes we have B^ cx y/C\D\ B^q (see Eq. (18) 
in Kleeorin et al., 2002). Note that for galaxies reasonable 
estimates are \D\ - 10 - 30, |C| ~ 0.1 and k - 0.3 - 0.5. 
With these parameters the present model gives toroidal 
field strengths of about the equipartition value. 

Fig. 2c demonstrates the change of sign of the effec- 
tive drift velocity with z, so that in one part of the disc 
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Fig. 2. As Fig. 1. The functions: a) rjAiz), b) c) 
Va{z). 



it is diamagnetic, and in the other it is paramagnetic. A 
diamagnetic velocity implies that the field is pushed out 
from regions with stronger mean magnetic field, while a 
paramagnetic velocity causes the magnetic field to be con- 
centrated in regions with stronger field. 

Fig. 6 shows the solutions for different values of the tur- 
bulent diffusion k of the magnetic helicity. It is apparent 
from Fig. 6 that the magnitude of the saturated toroidal 
magnetic field increases with k. The reason is clear; the 
increase of this parameter increases the flux of the mag- 
netic helicity, and causes a decrease of the magnetic part 
of the a-effect, thus increasing the total a-effect. 

When comparing the numerical and asymptotic solu- 
tions we need to take into account that the asymptotic 
solution (]2^)-(p^) was obtained only for k — 0. Thus such 
a comparison can only be performed for very small val- 
ues of K. If we compare the field at z = for, e.g., 
C — 0.1 and k ~ 0.05 we find that the difference between 
the asymptotic and numerical solutions is about 15 per- 
cent when D = —100. The novel feature, the maximum of 
the toroidal field i?^ at 2; > rather than z ~ 0, appears 
in the numerical solutions for k > 0.1, and is not described 
in the framework of the above asymptotic analysis. Note 



Fig. 3. As Fig. 1, with C = 0.1. 



also that a discrepancy between the numerical and asymp- 
totic solutions is perhaps not so surprising even for large 
values of \D\, as relation ( p4[ ) diverges near z — 1. 

In order to separate the study of the algebraic and 
dynamic nonlinearities we define a function Dy{B)/D = 
4>^{B) /[riA{B)riB{B)], using only the hydrodynamic part 
of the a effect. Thus, for B ^ l/-\/8 the function 
Dv{B)/D is a nonzero constant, because then rjA{B) ^ 
2/5/3, riB{B) - 3/5/3, ~ 2//32, with /3 = VSB. 

The saturation of the growth of the mean magnetic field 
in the case with only an algebraic nonlinearity present 
can be achieved when the derivative of the function 
dDv{B) / dB < 0. Thus, taking into account quenching of 
the turbulent magnetic diffusion we find that the algebraic 
nonlinearity alone (i.e. quenching of both the a effect and 
turbulent magnetic diffusion) cannot saturate the growth 
of the mean magnetic field because dDv{B) / dB for 
B > 1/V8. 

We will show here that the combined effect of the 
algebraic and dynamic nonlinearities limits the growth 
of the mean magnetic field. The dynamic nonlinear- 
ity is determined by the evolutionary equation (18) 
for x^- We introduce the parameter D]y{B)/D = 
a{B)/[a{B = 0)riA{B)r]B{B)] which characterizes both 
the algebraic and dynamic nonlinearities, while the pa- 
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Fig. 4. As Fig. 2, with C = 0.1. 

rameter Dv{B)/D characterizes only the algebraic non- 
linearity. The saturation of the growth of the mean mag- 
netic field is achieved when the derivative of the nonlin- 
ear dynamo number satisfies dDN{B) /dB < 0. We see 
in Fig. 7 that this condition is satisfied. However, we see 
also that dDv{B)/dB > 0. This implies that the algebraic 
nonlinearity alone (i.e. quenching of both the a effect and 
turbulent magnetic diffusion) cannot saturate the growth 
of the mean magnetic field. The same follows from the 
above asymptotic analysis. 

5. Simple models with radial extent 

Detailed numerical modelling of galactic dynamos is quite 
a complicated numerical problem. Dynamo models with 
conventional a-quenching are however very robust and al- 
low drastic simplifications that nevertheless reproduce ad- 
equately the basic features of galactic magnetic field struc- 
ture as reflected in the observational data. The aim of this 
section is to discuss to what extent these simplified models 
are comparable with the quenching discussed above. 

The basic no-z dynamo model for disc galaxies has 
proved to be a useful tool for studying dynamo action 
in these objects, and is described in Moss (1995). Here 



Fig. 5. The z-dependence of solutions for the local model; 
h/l = 5; C = 0.1, D = -40 and |x''(z)| = sm^{TTz/2). 
The various curves indicate results with k — 0.05 (dashed- 
dotted), K = 0.1 (dashed) and k = 0.3 (solid), a) the 
toroidal magnetic field B^{z), b) the poloidal magnetic 
field Br{z), c) the function x^(^)- 

we also include the tuning suggested by Phillips (2001), 
namely the multiplication by factors / = 7r^/4 of the terms 
representing the z-diffusion of B^. and B^. The no-z model 
differs from the local model of Sect. 2.1 in that it describes 



magnetic fields over the entire radial range, < ?' < i?, 
but all explicit dependence on the vertical coordinate z 
has been removed, with the first-order z-derivatives be- 
ing replaced by 1/h and the second-order z-derivatives 
being replaced by —f/h?. The field components Sr, B^ 
appearing in the no-z equations can either be thought of as 
representing mid-plane values, or as some sort of vertical 
average of values through the disc. 

For the no-z model in the axisymmetric case the mean- 
field dynamo equations take the form 



dB^ 
dt 



or 
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Fig. 6. The z-dependence of solutions for the local model; 
h/l = 5; C = -0.1, D = -40 and |x"(z)| = sin^(TTz/2). 
The various curves indicate results without K-quenching 
for K = 0.1 (dashed) and k = 0.3 (solid) and with k- 
quenching for k = 0.3 (dashed-dotted). a) the toroidal 
magnetic field B^(z), b) the poloidal magnetic field Bj.{z), 
c) the function x'^iz). 



Fig. 7. The dependencies of Dn{B)/D 
a{B)/[a{B = 0)r]A{,B)r)B{B)\ (solid) and Dv{B)/D = 
(f)v{B)l[rjA{B)'qB{B)\ (dashed) on the mean magnetic 
field for the one-dimensional model; h/l = 5; C = —0.1; 
K = 0.3; with different D: &) D = -20; h) D = -40; c) 
D = -100. 



dBr 



dt 



or 



dr 
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dr 
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where 



Va = [r]A 



s = - 



and A 



J__5 

j,2 Qj, 



{rBr) T^{rB^) + B^ 



d (I d 



dr \r dr 



{rBr) 



h/ R is the aspect ratio, Ra 
I Vr p\/ p and 770 is the maximum (unquenched) value of rj. 
The quenching functions (f>v , 4'm , VA, Vb and ryy contain in 
their arguments the factor 1 / ^/p because they are based on 
local equipartition at radius r. To determine the magnetic 
field distribution along the radius we use a Brandt rotation 
law, n{r) = no/[l + (r/r^)2]i/2 ^^^^^ ^ o.2, and the 



radial density profile p{r) = exp[~{r/rp)^ 

= 1. 



with Tp = 0.5, 

so that Ap — 2r/rp. We also set x"('') 

When comparing the no-z model with the local thin- 
disc model studied in the earlier parts of this paper, note 
that in the local thin-disc model, R^ = rh'^rj^^dfl/dr < 0. 
By the nature of the model, rdVL/dr is the value at a chosen 
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radius in the disc, and r does not further occur explicitly 
in the analysis. However the no-z model is global with 
respect to radius, and the value of rdQ/dr varies through 
the disc, from zero at r = to some maximum absolute 
value; for the Brandt rotation law this value is O.SIOq at 
r = ■y/2 / Svi^ . For the no-z model the global definition is 
= -Ra;(no — z) = ^oh^/ilT > 0. Less importantly, there 
are also small differences, of order 25%, in the effective 
values of Ra occurring in the two approximations, even 
though the formal definitions are the same - see Phillips 
(2001). 

We investigated solutions to these no-z equations for 
a range of parameter values. However we were unable to 
obtain satisfactory convergence to finite solutions without 
including a contribution to the diffusivity that was not 
quenched. Our feeling is that the simplest form of the no-z 
formalism may not be sufficiently robust to allow inclusion 
of the full 77-qucnching formalism described above. This is 
in sharp contrast to the situation without the inclusion 
of 77-quenching (Kleeorin et al. 2002), when satisfactory 
solutions with approximately equipartition strength fields 
were readily obtained. 

However, given the evidently satisfactory behavior of 
the local model presented in Sect. 4, it is apparent that 
extension to the radial range < r < i? in the manner 
described by Ruzmaikin et al. (1988, Ch. 7) would present 
no difficulties in principle, nor new features. 

6. Discussion and conclusions 

In this paper we present a more detailed description of a 
nonlinear galactic dynamo, that includes quenching of the 
turbulent diffusivity of the magnetic field in addition to 
the effects considered in our previous paper (Kleeorin et 
al. 2002). We find that as far as the model leads to results 
that are comparable with observations, these results are 
similar to those obtained from conventional galactic dy- 
namo models, with large-scale magnetic fields typically of 
equipartition strength, and with plausible values of the 
pitch angles. We confirm the conclusion of Kleeorin et 
al. (2002) that from a pragmatic point of view conven- 
tional models of nonlinear galactic dynamos are quite ad- 
equate to reproduce the directly observable manifestations 
of galactic dynamo action. 

Our approach is based on first principles as far as is 
possible in the framework of mean-field dynamo theory, 
and results in the conclusion that the self-consistent form 
of dynamo saturation is much more complicated than is 
suggested in conventional models for a galactic dynamo. 
We demonstrate the important role of two types of nonlin- 
earity (algebraic and dynamic) in the mean-field galactic 
dynamo. The algebraic nonlinearity is determined by a 
nonlinear dependence of the mean electromotive force on 
the mean magnetic field. The dynamic nonlinearity is de- 
termined by a differential equation for the magnetic part 
of the a-effect. This equation is a consequence of the con- 
servation of the total magnetic helicity. We have taken into 
account the algebraic quenching of both the a effect and 



the turbulent magnetic diffusion, and also dynamic non- 
linearities (see e.g. Rogachevskii & Kleeorin 1999, 2000, 
2001). Since the quenching of the a effect and the tur- 
bulent magnetic diffusion have the same origin (i.e. the 
direct and indirect modification of the electromotive force 
by the mean magnetic field), they cannot in general be 
taken into account separately. This implies there is no rea- 
son to include a quenching and to ignore the quenching 
of the turbulent magnetic diffusion, or vice versa. 

We have also verified that the algebraic nonlinearity 
alone (i.e. quenching of both the a effect and turbulent 
magnetic diffusion) cannot saturate the growth of the 
mean magnetic field. The situation changes when the dy- 
namic nonlinearity is taken into account. The crucial point 
is that the dynamic equation for the magnetic part of the 
a-effect (i.e. the dynamic nonlinearity) includes the flux 
of magnetic helicity. Without this fiux, the total magnetic 
helicity is conserved locally and the strength of the satu- 
rated mean magnetic field is very small compared to the 
equipartition strength. The inclusion of a magnetic he- 
licity flux means that the total magnetic helicity is not 
conserved locally because the magnetic helicity of small- 
scale magnetic fluctuations is redistributed by the helicity 
flux. In this case an integral of the total magnetic helic- 
ity over the disc is conserved. The equilibrium state is 
given by a balance between magnetic helicity production 
and magnetic helicity transport (see Kleeorin et al. 2002). 
Thus, the combined effect of algebraic and dynamic non- 
linearities limit the growth of the mean magnetic field and 
results in an equilibrium strength of the mean magnetic 
field which is of order that of the equipartition field, in 
agreement with observations of galactic magnetic fields. 

We find, perhaps quite naturally, that when includ- 
ing new physically significant effects we obtain less robust 
models. The limitations of simulations with such mod- 
els are connected not only with purely numerical prob- 
lems, which are however more severe than for conven- 
tional galactic dynamo simulations. In order to reproduce 
the detailed evolution of helicity and turbulent diffusivity 
there is the necessity for a deeper description of the real 
multiphase structure of the interstellar medium, and the 
physical processes that result in the development of heli- 
cal interstellar turbulence, than the standard description 
used in present-day models. Of course, the correspond- 
ing development of mean-field dynamo theory required to 
achieve this goal is far beyond the scope of this paper. 
We stress that direct numerical simulations of interstel- 
lar turbulence (see e.g. Korpi et al. 1999), followed by 
estimations of the appropriate mean-fleld dynamo control 
parameters, are very desirable and even essential in this 
context. 
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